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This work presents a reduced xt test to search for non-gaussian signals in the CMBR TT power 
spectrum of recent CMBR data, WMAP, ACBAR and CBI data sets, assuming a mixed density field 
including adiabatic and isocurvature fluctuations. We assume a skew positive mixed model with 
adiabatic inflation perturbations plus additional isocurvature perturbations possibly produced by 
topological defects. The joint probability distribution used in this context is a weighted combination 
of Gaussian and non-Gaussian random fields. Results from simulations of CMBR temperature for 
the mixed field show a distinct signature in CMB power spectrum for very small deviations (~ 
0.f%) from a pure Gaussian field, and can be used as a direct test for the nature of primordial 
fluctuations. A reduced xt test applied on the most recent CMBR observations reveals that an 
isocurvature fluctuations field is not ruled out and indeed permits a very good description for a 
flat geometry A-CDM universe, X930~ 1-5, rather than the simple inflationary standard model with 
X93o~ 2.3. This result may looks is particular discrepant with the reduced \ 2 °f 1-07 obtained with 
the same model in Spergel et al. (2003) for temperature only, however, our work is restricted to a 
region of the parameter space that does not include the best fit model for TT only of Spergel et al. 
(2003). 
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I. INTRODUCTION 

The new generation of cosmic microwave background 
radiation (CMBR) experiments has opened a new era 
in astrophysics, the precision cosmology era. Recent 
observations, especially those of WMAP (Wilkinson 
Anisotropy Microwave Probe), ACBAR (Arcminute Cos- 
mology Bolometer Array Receiver) and CBI (Cosmic 
Background Imager), brought a new light to CMBR fluc- 
tuations studies in large, intermediate and small scales. 
This new generation of experiments made possible the 
characterization of the power spectrum of temperature 
fluctuations up to the third acoustic peak (Hinshwaw et 
al. 2003; Pearson et al. 2002; Kuo et al. 2004). In- 
deed, the combination of CMBR and large scale structure 
(LSS) data allows cosmologists to constrain cosmological 
parameters for a given set of scenarios and also determine 
the nature of the primordial fluctuations. Since CMBR 
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carries the intrinsic statistical properties of cosmological 
perturbations, it is considered the most powerful tool to 
investigate the nature of cosmic structure. 

The most accepted model for structure formation as- 
sumes initial quantum fluctuations created during infla- 
tion and amplified by gravitational effects. The standard 
inflation model predicts an adiabatic uncorrelated ran- 
dom field with a nearly flat, scale- invariant spectrum on 
scales larger than ~ 1-2° (Guth 1981; Salopek, Bond & 
Bardeen 1981; Bardeen, Steinhardt & Turner 1983). Sim- 
ple inflationary models also predict the random field fol- 
lows a nearly-Gaussian distribution, where only small de- 
viations from Gaussianity are allowed (e.g. Gangui et al. 
1994). In hybrid inflation models (Battye & Weller 1998; 
Battye, Magueijo & Weller 1999), structure is formed 
by a linear combination of (inflation-produced) adiabatic 
and (topological defects induced) isocurvature density 
fluctuations. The topological defects are assumed to ap- 
pear during the phase transition that marks the end of 
the inflationary epoch. 

Correlated mixed field with adiabatic and isocurva- 
ture fluctuations has been considered by many authors, 
as Bucher et al. (2000), Gordon (2001) and Amendola 
et al (2002). Andrade et al. (2004) have suggested a 
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new mixed scenario describing a correlated mixture of 
two fields, one of them being a dominant Gaussian adi- 
abatic process to which is added a small contribution 
of an isocurvature field. In a previous work, Ribeiro, 
Wuensche & Letelier (2000) used such mixed scenario 
to probe the galaxy cluster abundance evolution in the 
universe and found that even a very small level of non- 
Gaussianity in the mixed field may introduce significant 
changes in the cluster abundance rate. Andrade, Wuen- 
sche and Ribeiro (2004) showed the effects of mixed mod- 
els in the CMBR power spectrum, combining a Gaussian 
(adiabatic) field with a second isocurvature field to pro- 
duce a positive skewness mixed density fluctuation field. 
This approach adopted a scale dependent mixture pa- 
rameter and a power-law initial spectrum to simulate the 
CMBR temperature and polarization power spectra for a 
flat A-CDM model, generating quite a large grid of cos- 
mological parameters combination. The results show how 
the shape and amplitude of the fluctuations in CMBR 
depend upon such mixed fields and how we can easily 
distinguish a standard adiabatic Gaussian field from a 
mixed non-Gaussian one and easily quantify the contri- 
bution of the second component. 

In this work, we apply a statistical test to both the 
mixed power spectrum simulations and recent CMBR 
data obtained by WMAP (Spergel et al. 2003), ACBAR 
(Kuo ct al. 2004) and CBI (Pearson et al. 2002) in or- 
der to compare how well the standard and mixed models 
describe the most recent observations. We also estimate 
the possible contribution of an isocurvature field to the 
primordial density fluctuation field in the mixed model 
scenario. 



logical parameter combination: Sltoth 2 = 1.02 ± 0.002; 
h = 0.71±°;^; fl C DMh 2 = 0.0224 ± 0.0009; n m h 2 = 
0.135±g;gg|; n^OmMpc- 1 ) = 0.93 ± 0.03 and a slope 
running spectral index with dn s /dln(k) — — 0.0311q'q]J 
(Spergel et al. 2003). The present CMBR data provides 
a strong support for adiabatic field domination. 

However, the WMAP data show several peculiarities 
at various values of 1. WMAP observations show that 
the fluctuations in large scale present a lower amplitude 
than the standard inflation model predicts. The tem- 
perature power spectrum is almost 30% suppressed on ~ 
1 degree angular scales (Hinshaw et al., 2003), specially 
those scales related to the quadrupolc (1 = 2) and the oc- 
topole (1 = 3), when compared with the predictions of the 
standard gaussian A-CDM models. Another impressive 
WMAP conclusions is the evidence of an optical depth 
to re-scattering by electrons with a value of r = 0.17 ± 
0.04, describing a reionization scenario at redshift 11 < 
z < 30 (Kogut ct al. 2003). 

In this paper, we try to reproduce the main feature 
of the CMBR power spectrum, in large, intermediate 
and small scales, invoking just the primordial physics, 
by considering a correlated mixed field in combination 
with six cosmological parameters: fib, fl c dm, ^A, n, H , 
the amplitude of the power spectrum, A, and the mixed 
coefficient, ao. Since we are interested in just cosmo- 
logical principles, no secondary effects are considered in 
CMB fluctuations, like reionization models, gravitational 
lenses nor Sunyayev-Zeldovich effects. 

III. THE MIXED MODEL 



II. RECENT CMB OBSERVATIONS 

The experiments ACBAR and CBI present the high- 
est sensitivity and highest signal to noise observations of 
CMBR temperature distribution in small angular scales 
( ~ 4 -5 arcminutes). In scales larger than the above 
mentioned (~ 20 arcminutes ) the WMAP satellite made 
a set of all-sky maps measuring both the temperature 
and polarization anisotropics, opening a new window for 
cosmological investigations. 

The WMAP observations, combined with LSS observa- 
tions has made possible to constrain a precise picture of 
the cosmos. Using the three-dimensional power spectrum 
P(k) from over 200,000 galaxies in the Sloan Digital Sky 
Survey (SDSS) in combination with WMAP and other 
data, Tegmark et al. (2003) show that recent observa- 
tions are consistent with a flat adiabatic A-CDM model. 
Specifically, the CMB power spectrum exhibits a first 
acoustic peak at 1 = 220.1 ± 0.8, with amplitude of the 
74.7 ± 0.5 [iK ; and a second acoustic peak at 1 = 546 ± 
10, with amplitude of 48.8 ± 0.9 fiK, a picture consistent 
with inflation predictions (Peiris et al. 2003). The best 
fit model compared with the WMAP plus ACBAR and 
CBI is for a A-CDM Universe with the following cosmo- 



A. The Mixed Primordial Power Spectrum 

In the mixed scenario, we suppose that the field has a 
probability density function of the form: 

V[5 k ] cx {l-a)h(5 k ) + af 2 {8 k ) (1) 

The first field will always be the Gaussian component and 
a possible effect of the second component is to modify 
the Gaussian field to have a positive tail. The parameter 
a in (1) allows us to modulate the contribution of each 
component to the resultant field. Like the hybrid infla- 
tion models (Battye & Weller 1998; Battye, Magueijo & 
Weller 1999), the mixture models consider the scenario 
in which structure is formed by both adiabatic density 
fluctuations produced during inflation and active isocur- 
vature perturbations created by cosmic defects during 
a phase transition which marks the end of inflationary 
epoch. 

Nevertheless, the mixed scenario considers a possible 
correlation between the adiabatic and the isocurvature 
fields on the post-inflation Universe. So, the fluctuations 
in super-horizon scales (inflated during the exponential 
expansion) are strictly uncorrelated and the mixing effect 
acts only inside the Hubble horizon, in sub-degree scales. 
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To allow for this condition and keep a continuous mixed 
field, a scale dependent mixture parameter, a = a(k) was 
defined. We assume the simplest choice of a(k), a linear 
function of k: 

a(k) = a k (2) 

To constrain the two component random field, we take 
6k = P(k)v 2 , where v is a random number with distribu- 
tion given by (f). So, we consider the primordial power 
spectrum of the mixed field in the form: 

P{k) mtx = M mix (a , k)P(k) (3) 

where the P(k) represents a primordial conventional 
power-law spectrum and M mlx (a) is the mixture term, a 
functional of fi and /g, which account for the statistics 
effect of a new component: 

M mx (a , k) ee f [{I - a k)h(v) + a kh(v)y dv (4) 

J V 

The effective mixed primordial power spectrum is ex- 
pressed as: 

P{k) mix =k n + M(a )k n+1 (5) 

where M(ao) represents only the coefficient dependence, 
cto- In the case of a pure Gaussian field, M(ao) ~ ao 
« 0, the mixed power spectrum will be a simple power- 
law spectrum, k n , where n is the conventional spectral 
index predicted by inflationary models. In the case of 
a mixed field, the phase correlations between both fields 
are estimated by the integral in Eq.(4), on mixture scales 
defined by Eq. (2). 



B. The Evolution of The Mixed Field 

To estimate the CMBR anisotropy we need to evalu- 
ate the evolution of fluctuations field generated in the 
early universe through the radiation-dominated era and 
recombination. In the mixed model, the evolution of both 
adiabatic and isocurvature components of the mixed den- 
sity field is considered an independent process, when only 
their effective amplitude correlation is considered on the 
decoupled surface of the CMBR photons. 

To compute the independent evolution of the adiabatic 
A-CDM mode and the isocurvature A-CDM mode, we 
have used the most precision code, the Linger function 
of the COSMICS code package (Bertschinger 1999; Ma 
& Bertschinger 1995). Linger does generate the most 
accurate results for the photon density field, being able to 
compute the CMBR anisotropy with integrations errors 
less than 0.15%. However, Linger routine is very time 
consuming, but for investigation of small deviations from 
gaussianity, the COSMICS package seems to be the most 
precise code for mixed CMB simulations. 



Once estimated the photon density evolution, the mul- 
tipole moments, Cj, of the CMBR temperature power 
spectra were estimated by a mixed photon density func- 
tion incorporated to the original COSMICS package: 



f k max 

Ci=4tt d 3 A;P™ a: (fc)(Ar ia; ) 2 (fc,T) (6) 

Jo 

The function A™ 2 represents the mixed photon density 
field in the last scattering surface defined by: 



A? lx = {1- a k)Af dl + a fcAf so (7) 

where Af dl and Af so are the photon density function 
estimated by COSMICS, respectively, for an adiabatic 
and an isocurvature seed initial conditions. 

Inserting (7) in (6), we obtain a mixed term in the Cj 
estimation. This condition suggests that the amplitude 
of both fields are cross-correlated at the last scattering 
surface, with a mixing ratio defined by «o and in a char- 
acteristic scale defined by a k. The power spectra esti- 
mated by Eq.(3) consider a flat A-CDM universe slightly 
distorted by a non-gaussian statistics with constant spec- 
tral index. 

A number of realizations of CMBR temperature power 
spectra were made and the mean temperature fluc- 
tuations combining Gaussian, Exponential, Lognormal, 
Rayleigh, Maxwcllian and a Chi-squared distributions 
was estimated (Andrade et al. 2004). The simulations 
show that the influence of the specific statistics of the 
second component in the mixed field is not so important 
as the cross-correlation between the amplitudes of both 
fields. However, some important results were obtained. 
It is possible to directly assess and quantify the mixture 
of a correlated adiabatic and isocurvature non-Gaussian 
field. This behavior points to another possibility to ex- 
tract information from a CMBR power spectrum: the 
possibility of detecting weakly mixed density fields, even 
if we can not exactly identify the mixture components 
distribution. 

In Figure 1, we see how the shape and amplitude of 
the spectrum change even for small values of a (~T0~ 4 - 
10~ 3 ). The peaks intensities are clearly susceptible to the 
existence of mixed fields, although distinguishing peaks 
of higher order (2 nd , 3 rd , etc.) in the mixed context is not 
a straightforward task, since their intensities, compared 
to the first peak, are very low. The effect of increasing ao 
is a power transfer to smaller scales (1 >I000), while the 
super-degree scales are less affected. However, for large 
values of ao (ao >3xI0 -3 ) there is a fast increase in the 
temperature fluctuations, probably caused by correlation 
excess between the mixed fields, resulting in more power 
in small scales. Therefore, an acceptable range for a is 
set to be a < 3 x 10~ 3 . 
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FIG. 1: CMBR temperature mixed angular power spectrum 
estimated for a A-CDM model in different mixing degrees. 

IV. A x 2 TEST ON CMB OBSERVATION 

In order to quantify a possible non-gaussian isocurva- 
ture contribution to the fluctuations field, we applied a 
maximum likelihood approach to the power spectrum es- 
timated by the most recent CMBR observations. We ran 
a set of CMBR realizations for mixed and pure density 
fields, considering a flat A-CDM Universe consistent with 
inflation predictions and LSS observations. 

Instead of the usual CMBR analysis approach, which 
considers tensors contribution and reionization effects, 
we consider only the main six cosmological parameters 
(ujcdm, Mb, u\, n, Ho, ao), and also A, described in the 
range of possibilities: 0.8< n <1.2; 0.015 <n b < 0.03; 0.6 
<Ha< 0.8; 60 <Ho< 80; cold dark matter density set as 1- 
(flb + Ha), ranging of 0.170 < Q c dm< 0.385; and the range 
values of ao se t in 0.0 < a$ < 3xl0~ 3 . With this con- 
siderations, we ensure the simulations are in agreement 
with the basic preferences of standard inflation scenar- 
ios: flat geometry (fltot ~ 1) and a nearly scale invariant 
primeval spectrum (n ~ 1). f2& =(0.019 ± 0.01)h~ 2 is 
consistent with the mass density of baryons determined 
by Big Bang nucleossynthesis and the large scale struc- 
ture observations, which suggests that the Hubble con- 
stant H assumes values in the range (66 < H < 75) 
km/sec. Mpc; and a positive cosmological constant value 
in the range (0.065< J7a <0.75), as estimated recently by 
Tegmark et al. (2003). The amplitude of the primordial 
power spectrum, A, was estimated by the best fit value 
for each model simulated. 

We applied a reduced \ 2 test (Bevington 1969) be- 
tween the CMB data and the power spectra simulations 
in the mixed context in four steps. First, we consider the 
WMAP, CBI and ACBAR data independently, estimate 
the reduced % 2 for all three data sets and then add all 



data to obtain a new estimation of the reduced x 2 - So 
we estimate the residual difference for C; estimation in 
^K 2 , and consider the power spectrum data points (899 
for WMAP; 24 for CBI and 14 for ACBAR) minus 7 de- 
gree of freedom. We create a grid of nearly one hundred 
simulated values of ao for each combination of standard 
cosmological parameters. Since COSMICS is very time 
consuming, the variations in the standard cosmological 
parameters were set in order to give just the main direc- 
tion for the x 2 face variations in cosmological parameters. 
The most precise models, as defined by the x 2 test, are 
summarized on Table 1, which also contains the best ao 
value estimated for each parameter combination listed 
with 68% confidence limit. 

The simulations clearly show that WMAP data tends 
to favor mixed models. The x 2 estimated for the fit be- 
tween mixed models and WMAP data are significatively 
lower for mixing degrees of order ~10 -3 , in all combi- 
nation of cosmological parameters tested. However, CBI 
and ACBAR data in small angular scales tends to favor 
the adiabatic standard models. The combination of all 
data are more sensible to WMAP influence (more signal 
in large and intermediate scales with low error bars) re- 
sulting in best fits models in a mixed scenario. All data 
favor a A-CDM Universe with contents of barions, dark 
matter and dark energy coherent with that estimated by 
the WMAP and LSS analysis. 

Figure 2 shows the powers spectrum of the standard 
adiabatic and mixed models for the three best fit models 
obtained by the \i test. It also shows the behavour of 
the all data xt while varying the ao value. It is clearly 
seen that the mixed model adjust the data better than 
the standard model in basic combination of cosmologi- 
cal parameters. In this figure, is possible to observe that 
the relation between the amplitudes of the acoustic peaks 
in the mixed scenario is more adequate to describe the 
CMBR power spectrum feature without consideration of 
any model of reionization or tensors contributions. Al- 
though the mixed model predicts higher quadrupole and 
octopole contributions that predicts the standard adi- 
abatic A-CDM model, the plots clearly show that the 
mixed predictions are consistente with acoustic peaks 
amplitude estimated by WMAP. Indeed, a model estima- 
tion in a good fit is obtained for a strictly flat Universe, 
rather than an open Universe or, instead, a slope running 
spectral index (Kosowsky & Turner, 1995) as suggested 
by the WMAP team. 

One possible explanation of the low amplitude in large 
scales and the high amplitude of the secondary acous- 
tic peaks predicted by the standard model (pure adia- 
batic) simulated in this work is the lack of the high opti- 
cal depth, as considered by the WMAP team (Hinshaw 
et al. 2003; Spergel et al. 2003). Reionization affects 
the CMBR in several manners, through Thomson scat- 
tering of photons from free electrons in the intergalac- 
tic medium, which suppress the primordial anisotropics 
of the CMBR imprinted on the decoupled surface. In 
a reionizing scenario, the acoustic peaks in the CMBR 
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"Best model estimated for WMAP data. 

6 Best model estimated for CBI data. 

c Best model estimated for ACBAR data. 

d Best model estimated for all data (plot on Figure 2A ). 

^Second model estimated for all data (plot on Figure 2C). 

^Third model estimated for all data (plot on Figure 2B). 

TABLE I: The reduced \ 2 test applied on the WMAP, CBI and ACBAR data and the power spectrum simulations in a mixed 
scenario for some A-CDM models. The last column shows the best ao estimation for each model with 68% confidence limit 



power spectrum would have the amplitude lowered by 
e~ 2r and the rescattering would generate large angle 
temperature fluctuations through the Doppler effect (Hu 
& White, 1997). In this work, rather than considering 
rcionization effects and invoking just primordial physics, 
we can explain the main feature of the CMBR power 
spectrum in large, intermediate and small scales in the 
context of mixed density field. 



V. DISCUSSION 

We applied a modified statistical xt test to the most 
recent CMBR observations. The results allow us to com- 
pare the performance of power spectra simulated with 
mixed and standard density fields to describe the mean 
feature of CMBR temperature fluctuations. The results 



showed that the WMAP data tends to favor mixed den- 
sity fields with isocurvature fluctuation contribution of 
order ~10 -3 (xig2^ 1-5 for mixed model against Xs92~ 
2.4 in standard models) without the need of considering 
any model of rcionization. Also, large scale fluctuations 
are in good agreement with WMAP measurements. Our 
result may looks discrepant with the reduced \ 2 of 1.07 
obtained with the same model in Spergel et al. (2003) 
for temperature only, however, our work is restricted to 
a region of the parameter space that does not include the 
best fit model for TT only of the WMAP team. We can 
conclude that the mixed scenario offers a good alterna- 
tive to describe the shape of TT CMBR power spectrum. 
Indeed, the data favors a A-CDM Universe with barions, 
dark matter and dark energy contents which are coher- 
ent with a strictly flat Universe, with no need to include 
tensor contributions nor a slope running spectral index. 
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FIG. 2: CMBR temperature power spectrum estimated for the best three parameter combination: flt= 0.03; Qcdm= 0.27; 
n A = 0.7. In A: H = 65Km/sMpc and n=l.l; in B: H = 70Km/sMpc and n=1.2; in C: H = 70Km/sMpc and n=l.l. The 
small picture in figures show the behavior of \t while varying ao value 
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For a A-CDM Universe with contents of flb= 0.03; and 
the best value indicated for the spectral index, by the x 2 
test, n~l.l; the contribution of the isocurvature field is 
estimated as a = 0.00042 ± 0.0003 with 68% confidence 
limit. 

Although the standard A-CDM model fits well a small 
number of data points of CBI (x 2 4~ 1 f° r standard mod- 
els against X24^ 2.3 in mixed models) and ACBAR (x?~ 
2.4 for standard models against xf~ 6.8 in mixed mod- 
els), when considering all data points, CBI, ACBAR and 
WMAP, the mixed model clearly shows ability to de- 
scribe the mean feature of TT CMBR observation in 
large, intermediate and small scales. With the upcoming 
generation of experiments which will probe temperature 
fluctuations in small scales, such as Planck satellite, it 
will be possible to set new, more stringent, estimates of 
the parameter spaces and isocurvature fluctuations con- 
tributions. Also, we will try to constrain the mixed model 
predictions upon the temperature-polarization cross cor- 
related power spectrum, in order to consider reinoniza- 
tions models in the mixed scenario. Presently, we con- 
clude that isocurvature fluctuations can not be ruled out 
in cosmological studies. Also, conflicting with the results 
of WMAP map team (Komatsu et al. 2003), the re- 
sults of some independent searches for non-gaussianity in 
WMAP maps evidence that WMAP data are consistent 
with gaussian condition for 1 <250, marginally gaussian 



for 224<1 <350, and non-gaussian for 1 >350, as revealed, 
for instance, by the phase mapping technique (Chiang et 
al. 2003; Coles et al. 2003). This picture points an extra 
expectation in the study of the mixed model predictions, 
altough, in order to constrain a better estimative of ao, 
rcionizations models must be investigated in the context 
of mixed fields. 

The next step in the investigation of the mixed mod- 
els in CMBR features will be the predictions of the 
reionization effects, even in a low redshifts, the degen- 
eraices between cosmological parameters in mixed con- 
text, the predictions of the angular correlation function 
and the temperature-polarization cross correlated power 
spectrum. 
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